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Abstract 

Previous computations of the potential landscape with the shapes parameterized 
in terms of Cassini ovaloids are extended to collective dynamics at hnite excitations. 
Taking fission as the most demanding example of large scale collective motion, 
transport coefficients are evaluated along a fission path. We concentrate on those for 
average motion, namely stiffness C, friction 7 and inertia M. Their expressions are 
formulated within a locally harmonic approximation and the help of linear response 
theory. Different approximations are examined and comparisons are made both with 
previous studies, which involved different descriptions of single particle dynamics, 
as well as with macroscopic models. Special attention is paid to an appropriate 
definition of the deformation of the nuclear density and its relation to that of the 
single particle potential. For temperatures above 2>MeV the inertia agrees with 
that of irrotational flow to less than a factor of two, but shows larger deviations 
below, in particular in its dependence on the shape. Also friction exhibits large 
fluctuations along the fission path for small excitations. They get smoothed out 
above 3 — AMeV where 7 attains values in the range of the wall formula. For 
T > 2MeV the inverse relaxation time f3 = 'j/M turns out to be rather insensitive 
to the shape and increases with T. 

PACS numbers: 21.60.Ev, 21.60.Cs, 24.10Pa, 24.75+i 

to appear in PRC 


1 Introduction 


One of the oldest but still most challenging problems of nuclear physics is an adequate 
description of collective motion at hnite excitations. As the prime example one may 
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quote nuclear fission which has attracted the attention of both experimentalists as well 
as theoreticians since its discovery. To date it is still an open question which type of 
conhgurations the system undergoes on its way from the potential minimum over the 
saddle region down to scission. Whereas in the early days those of the compound model 
were clearly favoured in theoretical pictures, after the discovery of the shell model that 
of independent particle motion came into fashion more and more. This development 
was enhanced after computers got fast enough such that Hartree-Fock type computations 
could be done in every lab. 

However, there can be little doubt that this picture fails to describe collective motion 
at hnite excitations where one is compelled almost by experimental evidence that the dy¬ 
namics shows irreversible behaviour, not only by the very nature of the decay process itself 
but by the appearance of frictional forces. It is more than questionable that this feature 
can adequately be met by introducing simple minded collision terms. Decent descriptions 
of hssion in terms of the one body density operator most likely require to consider cor¬ 
relations beyond the independent particle picture, together with non-Markovian effects. 
This is a difficult problem in itself, not to mention the computational task of solving this 
equation of motion for the one body density. 

For these reasons it may still be interesting and worth while to start from a more 
phenomenological point of view introducing the shape parameters as collective variables. 
It is true that in this way again the picture of independent particles will serve as a 
starting point, in the form of the deformed shell model. However, the latter is simple 
and flexible enough to allow one considering residual interactions, in one way or other. 
As we shall see, it may be possible to gain insight into their importance by studying 
dynamical aspects. Likewise, we may be able to get information on the complexity of 
the conhgurations which are to be considered. Such a task becomes more feasible in case 
collective motion is sufficiently slow. Then one may exploit the quasi-static picture which 
reduces the complexity of the full problem drastically. Under such circumstances one 
may actually linearize the problem and treat collective motion locally within a harmonic 
approximation. In this way one may take advantage of the benehts of linear response 
theory. 

One of the major problems in theories of this type is to hnd a decent guess for the 
relevant macroscopic variables, a problem which is familiar almost from all transport 
theories. For nuclear hssion there exists some kind of guiding principle through the liquid 
drop model. The latter is known to represent the static energy for temperatures above 
1 — 2 MeV. Since at these temperatures one expects motion to be strongly damped, 
it will most likely follow somehow the line of steepest decent. Possible shapes which 
a hssioning nucleus may assume on its way to scission have been looked for in |l[] by 
minimizing the liquid drop energy. This minimization has been done for some realistic 
energy density functional under the constraint of hxing a parameter which measures the 
distance between the evolving fragments. Incidentally, it is the same parameter which we 
are going to exploit later on in our approach. It so turns out that the shapes found in this 
way can be approximated fairly well by the Cassini ovaloids introduced to nuclear physics 
in [Q]. Later in |Q] a single particle model has been constructed for such a parameterization 
of shapes, which was based on the Woods-Saxon potential. 

In this paper we are going to use this model for computations of transport coefficients, 
after some suitable modihcations which are necessary to incorporate the effects mentioned 
above. One of our goals will be to study average motion along the fission path for different 
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temperatures, as it is reflected in the associated transport coefficients of inertia, friction 
and local stiffness. In this sense the aim of our present work is similar to the one of 0, 
where a two center shell model was used. The latter feature renders the previous model 
simpler on the computational level. On the other hand, the parameterization of the shape 
by means of Cassini ovaloids opens the possibility of treating more realistic shapes, which 
are perhaps better suited to describe the later stages of a hssion process. Furthermore, 
it is fair to say that the Woods-Saxon potential may be supposed to resemble more the 
’’true” mean held. 

For Cassini ovaloids commonly a few parameters suffice to treat in simple terms a whole 
variety of realistic shapes including very compact ones as well as strongly deformed ones 
with a well developed neck, or even those corresponding to separated fragments. In this 
sense this parameterization may be considered superior to expansions in terms of spherical 
harmonics (see 0). In the ideal case one would then be able to compute transport tensors 
for all the parameters, the collective degrees of freedom, one claims to be relevant. This is 
a tremendous task and so far has been carried through only for a two dimensional model 
[^, without utilizing though the full microscopic potential of linear response theory. In 
this paper we want to restrict ourselves to the one dimensional case. The main reason for 
that is found, of course, in the simplihcation one gains by this restriction. However, it 
may be said that at present most of the applications of macroscopic equations of motion 
to hssion at hnite excitation adhere to a similar conhnement, see e.g. 0-[0- Evidently 
one then needs to rely on the ’’right” guess of the hssion path. As said before and for 
arguments given there we presume it to be represented well enough by the line along the 
valley of the static energy. Possible improvements have to be left for future studies. 


2 Deformed shell model 


2.1 Shape parameterization 

We follow the suggestion put forward in , but would like to repeat the most important 
elements for convenience. The Cassini ovaloids are obtained by rotating the curve 


p{z, e) = Ro 


+ Aez'^/R.Q 



( 2 . 1 ) 


around the z-axis, with z and p being cylindrical coordinates. The constant a is dehned 
by volume conservation, implying that the family of shapes (|2.1| ) depends only on one 
deformation parameter e. As is easily recognized from (|2.1|) the value of e = 0 corresponds 
to a sphere. For 0 < e < 0.4 the form resembles very much that of a spheroid with the 
ratio of the axes given by 

shorter axes 1 —2e/3 

--=- 2.2 

longer axes l + e/3 

At e ~ 0.5 a neck appears and at e = 1.0 the nucleus separates into two fragments. A few 
examples of the family (|2.1|) are shown in Fig.l. 

It is possible to describe a more general class of axially symmetric shapes by exploiting 
an expansion about the surface given by ( |2.1| ). We may introduce two new coordinates R 
and X such that R = const = Rq corresponds to the Cassini ovals (|2.1|) . The coordinate x 
specihes the position of a point on the line given by ( 0 . see for details. With these 
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two variables at our disposal we may parameterize a new shape. The latter is meant to 
express the deviation from the ovaloid given by ( p.l|) by means of an expansion into a 
series of Legendre polynomials 

R{x) = Ro{l + J2 anPn{x)) (2.3) 

n 

The full set of collective variables or parameters then includes the coefficients in addi¬ 
tion to the e from before. 



Figure 1: The parameterization of the shape of the nuclear surface in terms of Cassini 
ovaloids. The values of the deformation parameter e (see (|2.1|)) are indicated in the Figure. 


Sometimes it is convenient to introduce a measure for the overall elongation of the 
nucleus instead of e. One may chose, for instance, the distance R 12 between the left 
and right center of masses. To have a dimensionless quantity one may divide R 12 by the 
diameter 2Rq of the sphere (of identical volume) to get 


_ Ri2 I ^cm I dV 

“ m ;; “ R„sdv 


(2.4) 


with Zcm being the ^-coordinate of the center of mass of the the whole complex. The 
integration is carried out over the volume within the sharp surface specihed by (2.3). 
Asymptotically, the ri 2 turns into half of the distance between centers of mass of the 
hssion fragments. Incidentally, this variable is dehned uniquely for any parameterization 
of the shape and has been used in the past by many authors. This feature facilitates 
comparisons to theories which are based on shell models with different shape variables. 
The relation between ri 2 and e is demonstrated in Fig.2. 


2.2 Single-particle Hamiltonian 

The single-particle Hamiltonian hipm will be constructed like in [^. It has terms for the 
kinetic energy T, the radial potential V, the spin-orbit coupling I 4 . 0 . and the Coulomb 
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potential Vcoui- 


(2.5) 


hipm — T + V + Vs,o. + Vcoul 


The radial part V is represented by a finite depth Woods-Saxon potential 

V{p,z) = Vo[l +exp{l{p,z)/a)]-^ 


( 2 . 6 ) 


where l{p, z) is the shortest distance from the point (p, z) to the sharp surface and a is the 
diffuseness parameter which is assumed to be constant along the surface. The spin-orbit 
potential may be written in a way which makes apparent that it is proportional to the 


gradient of the potential given in (|2.6|) , 

Vs.o. oc [s,p\VV 


(2.7) 


Here p and s stand for the nucleon’s momentum and the spin. The Coulomb potential is 
calculated for a charge distributed uniformly inside the sharp surface (p.3|) or (H). 

The single-particle energies and wave functions are determined by diagonalizing the 
matrix of the Hamiltonian (|2.5|) calculated with the wave functions of a deformed axially 
symmetric oscillator potential, see m- An example of the deformation dependence of 
the single-particle energies is shown in Fig.3. As the result of diagonalization one obtains 
not only energies and wave functions of bound states but also those of discrete states of 
positive energy, which for the Woods-Saxon potential lie in the continuum. The density 
of these states depends on the number of oscillator shells included into the basis. In 
the computations within the shell correction method the number of oscillator shells is 
optimized by the requirement that the states with positive energies provide a smooth 
extrapolation of the density of bound states into the continuum. Accounting for such 
states with positive energy improves considerably the ’’plateau” of the shell correction 
as a function of the averaging interval. In the present paper we do so not only when 
calculating deformation energies but also in the computation of transport coefficients. 
The cut-off in the single-particle energy was set equal to 2QMeV. We have checked that 
a variation of the cut-off energy within the interval 5 — 20MeV does not change much 
the values of the transport coefficients. This may be understood from the fact that these 
transport coefficients refiect the truly low frequency behavior of the system. 


2.3 Deformation energy 

The deformation energy Edef at zero temperature is calculated according to the shell 
correction method |]T3,[0 as the sum of the liquid drop energy and the shell 

correction SE'^’P + (including the one for the pairing energy) 


Edef — E. 


LDM 

def 


J2{SEP’^ + 6PP’^) 


( 2 . 8 ) 


p,n 


The liquid drop energy is computed as the contributions from the Coulomb EqouI and 
surface Es energies according to ||l2[|,[p!^ 


rpLDM _ \ TP 

^def - ^Coul + 


+ EPs) 


(2.9) 


where Eq^^i and Eg are the corresponding energies of the spherical shape. As an example. 
Fig.4 exhibits the results of the calculation of Edef and E^^^ at zero temperature for 
as function of parameters e and a^. 
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The temperature dependence of Coulomb and surface energy is accounted for by using 
the forms 


Ecoui{T) = EcouiiT = 0)(1 - «T2), Es{T) = Es{T = 0)(1 - (3T^) (2.10) 


with a = 0.000763MeC ^ and {3 = 0.00553 MeV ^ [^. To compute the shell correction 
at hnite temperature we use the phenomenological ansatz proposed in [ffl 


SEP'^{T) + 6PP'^{T) = [5EP'^{T = 0) + 6PP'^{T = 0)] 


r 


sinh T 


( 2 . 11 ) 


with T = 2Ti'^T/hujQ and hcuo = A For temperatures larger than IMeV (as 

considered in the present paper) the shell effects are strongly suppressed. Already at 
T = IMeV the minimum of the total deformation energy almost coincides with the 
bottom of the liquid drop valley. As said before the latter can be approximated rather 
well just by the Cassini ovaloids. Thus we restricted our set of deformation parameters 
to the one parameter e only, with all the put equal to zero. However instead of e we 
prefer to use dehned by |2.4| . The ri 2 -dependence of the total deformation energy and 
that of the liquid drop is shown in Fig.5. being the mass number of the hssioning nucleus. 


3 Dynamics in locally harmonic approximation 


In the following we will assume to be given a many body Hamiltonian H{xi,pi,Q{t)) 
which depends parametrically on the collective variable Q which specihes the shape of 
nuclear surface. Although for the computations to be presented below, it will mostly be 
identical to ri 2 , in this section we still prefer to use the general notation Q instead, last 
but not least to indicate the general validity of the discussion to come. This Hamiltonian 
is assumed to represent the system’s total energy. On the level of the shell model this 
means to add some c-number terms to the sum over those single particle Hamiltonians 
introduced in ( |2.5| ) (see[]^). As will be discussed in the next section, later on we want 
to account for collisional damping, which from a principal point of view requires to add a 
two body interaction Vj3}{xi,pi). For the moment it is not very important to know details 
about the way it will be handled — besides the fact that we claim this interaction to be 
independent of Q. 

As a consequence of the latter feature, the generator for collective motion, namely 
dH{xi,pi,Q)/ dQ = F{xi,pi,Q), is of pure one body nature. This operator dehnes 
the main source of the coupling between the collective degree of freedom Q(t) and the 
nucleonic ones. Indeed, within the locally harmonic approximation (LHA) the effective 
Hamiltonian can be written as 


H{Q) = H{Qo) + (g 


Qo)E + -(g 


go)^ 



(3.1) 


QoTo 


In the second order term the ’’nucleonic” part appears only as an average of the corre¬ 
sponding operator. Consistently with the harmonic approximation, this average is to be 
built with that density operator Pqs(go) which in the quasi-static picture is to be cal¬ 
culated with the Hamiltonian at Qo, namely H{Qo)- It is here where thermal concepts 
come into play. In this ’’unperturbed” density operator (for the nucleons) one needs to 
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specify the amount of heat the nucleonic (or intrinsic) degrees of freedom have at the given 
conhguration parameterized by Qq. The simplest possibility is offered by the canonical 
distribution []Pqs(Qo) cx: exp(—II(Qo)/T), to which our computations will be restricted. 
Clearly this picture has to rely on the assumption of a quick relaxation of the relevant 
internal degrees of freedom; we will come back to this question later on. 

0,113 


Details about this LHA can be found in many references, see e.g. 


|T8| : There 


it is also described how this local dynamics can be handled within a suitable application of 
linear response theory. For this reason, we will only recall the most important theoretical 


issues. 


3.1 Collective response function 

The local motion in the Q-variable can be described in terms of the so called collective 
response function Xcoii(<^)- It can be derived by introducing a (hypothetical) external 
force Ffext{t) and by evaluating how the deviation of {F)^ from some properly chosen 
static value reacts to this external held in linear order: 


Xcoll{uj) f 


(3.2) 


As shown in and the Xcoii(a.') can be brought to the form 


x(^^) 

Xcoll(ij) ~ 1 I 7 / \ 

1 + kx{uj) 


(3.3) 


Here a response function x(a;) for ’’intrinsic” motion appears: The x(co’) measures how, 
at some given shape Qq and for some temperature Tq, the nucleonic degrees of freedom 
react to the coupling F5Q{ui). Its time dependent version reads 


Xit - s) = Q{t -s)^ tr(pq^(Qo, To) [F(t), F(s)]) = 2iQ{t - s)x!'{t - s) (3.4) 

In this expression the time development of the held operators is dehned by the same 
Hamiltonian H{xi,Pi,Qo) which appears in the density p^s- The function x"(t — s) on 
the very right stands for the so called dissipative part. In Fourier space the full response 
separates into real and imaginary parts like x(a;) = x'(a;) + ix"(a;), the x'(^) sometimes 
being called reactive part. 

As one may guess from the very construction, the derivation of (|3.3|) relies on quasi¬ 
static properties of the nucleonic degrees of freedom. For instance, such properties appear 
in the coupling constant k which is to be determined from 


- k-^ 



0) = C(0 )+x(0) 


(3.5) 


Moreover, the nucleonic degrees of freedom are assumed to behave ergodic in the sense 
of having the adiabatic susceptibility x^"^ = ~^{F)/5Qs be identical to the isolated one, 
the static response x(0): 

x(0) = (3.6) 


^For a discussion of the general problems of using the concept of temperature for an isolated system 


see I IF 
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As has been demonstrated in this condition is not fulhlled in the deformed shell model, 
which implies that special measures are to be taken to which we will come to below. 

It should be noted that the derivation of (|3.3|) involves a self-consistency relation 
between the deformation Q of the mean held and the one of the density. The latter may 
be measured by the expectation value For linearized dynamics this self-consistency 

condition reduces to the equation 




(3.7) 


well known from the case of undamped vibrations i, El- Realize, please, that the 
quantity {F)^ is to be calculated with the actual dynamical nuclear states accounting for 
their appropriate occupations. 

Before concluding this subsection we like to write down a more convenient form for 
the stiffness C'(O) appearing in ( p.SI) . It is dehned as the second derivative of the internal 
energy E{Q, S) with respect to deformation at hxed entropy S. Since it is not easy to 
calculate such a derivative it is better to reexpress it by the one of the free energy / at 
hxed temperature (see eqs.(A.18)-(A.19) of 0) 


C(0) = 


fd^E\ 

[Wj, 


'm + 

OQ^t 


aV(Q,T) 

^dQ) g dQdT 


(3.8) 


In [17] it was found that for temperatures larger than 1.5 — 2.QMeV the change of T 


with the collective coordinate Q is small such that the second term on the right hand 
side of eq. (|3.8|) can be neglected. Below we would like to use the shell correction method 
when calculating static energies. In this method usually the intrinsic energy is involved, 
rather than the free energy. Thus one needs to relate the derivatives of the free and the 
intrinsic energy both taken at at hxed temperature. This can be done by differentiating 
the relation E = f + TS with respect to deformation and obeying that the entropy can 
be expressed as S' = —{df/dT)Q. As the result one gets 


^ ’ \8Q^ L \8Q'^ L dT\aQ-‘ L \aQ\ 


(3.9) 


Like before, the expression on the very right is justihed for larger temperatures, where 
the change of d‘^f/dQ‘^ with T is small. This approximation is used in the computations 
presented below. In Fig.6 the stiffness (|3^) of the energy (^^) - (p.ll|) is shown as function 
of the deformation parameter ri 2 . It is seen that its liquid drop part becomes negative for 
ri 2 > 0.5 and the total stiffness exhibits rather strong huctuations due to shell structure. 


3.2 Transport coefficients 

In general the frequency dependence of Xcoii(<^) exhibits a complex structure. Its dissipa¬ 
tive part Xcoii(‘^) represents the strength distribution over all possible modes of the whole 
nucleus which can be excited by an external force like the one introduced above, namely 
Ffext{t)- Rewriting (|X^) as 




(3.10) 





















it follows that the inverse of the collective response function can be interpreted as an 
integral kernel for the effective equation of motion for the time dependent quantity S{F)t. 
Evidently, this time dependent form of the equation of motion must be expected to contain 
non-Markovian effects. However, there may be situations for which it becomes possible 
to reduce this complicated structure to differential form. The clue to this simplihcation 
can be found by recalling the case of the damped oscillator for which ( p.lOj ) takes on the 
form 

(Xosc(o;))“' 6{F)^ = + C^) 6{F)^ = -feA^) (3.11) 

Here the following effective ’’forces” appear: an inertial one, a friction force and the 
conservative one which is related to the derivative of the effective collective potential, 
calculated in linearized form. The associated ’’transport coefficients” have been marked 
by an upper index ” F”, to indicate that they are associated to the quantity 6 {F). 

Now the strategy of how to handle the general case is clear: Whenever there is a pro¬ 
nounced peak in the strength distribution AowA) may approximate it by a Lorentzian. 
Since the latter is dehned by three quantities, like the width, the position of the maxi¬ 
mum and its height, one may deduce the three transport coefficients which appear in the 
oscillator response, for instance as given by (|3.11|) . In principle, such a procedure may be 
applied to any one of the peaks of the original strength distribution, which allows one to 
deduce these transport coefficients for all the possible collective modes, the low frequency 
ones as well as the high frequency ones. As a matter of fact, such a scheme is commonly 
adopted for the description of collective vibrations, as discussed in . The difference 


to the present application is seen in the fact that we want to apply this procedure to de¬ 
scribe global dynamics in locally harmonic fashion. The latter aspect puts an additional 
constraint. As mentioned earlier, our application of linear response theory goes along with 
the assumption that collective motion is slow, such that the nucleonic degrees of freedom 
follow closely a thermal equilibrium. We may recall that in some formulas given above, 
this equilibrium has been parameterized by the quasi-static density operator pqs(Q,T). 
Estimates of the time scale on which such a relaxation can be expected to occur will be 
given below. But already on this level of information, it is clear that the whole concept 
would probably not work if high frequency collective modes would be important. For 
this reason we are bound to concentrate on low frequency ones, if the construction of the 
transport coefficients is to be consistent with the basic assumptions for the applicability 
of the quasi-static picture. 

In practice the ”£t” of the strength distributions of the oscillator model to the ’’correct” 
one involves the full response functions, not only their dissipative parts. In short such 
an adjustment may be characterized by (Xcoii(<^))~^ ^A)u) — (Xosc(‘^))~^ A more 

correct form can be written as a variational procedure 


k'^XcoiiA) 


—— i'yAiA + C{uJi) 


du = 0 


(3.12) 


with the variation to be performed with respect to the coefficients C{uJi), xAi) and 
M{uJi). In this way both real and imaginary parts of Xcoii(<^) are fitted simultaneously. 
Fortunately, for the practical applications to be discussed below it turns out that the 
coefficients ^(a;!), 7 (a;i) and M{uji) are rather insensitive to the upper integration limit 
^max- The value of the latter was hxed to be humax = 5MeH. 


The reader will have noticed the appearance of the factor in (|3.12|) , as well as the 


fact that we have left out the upper index F in the transport coefficients. This change is 
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easily understood by referring to the self-consistency condition ( |3.7|) . The latter can be in¬ 
terpreted as a transformation from the F-mode to the Q-mode. It implies a corresponding 
transformation both for the response functions as well as for the transport coefficients (for 
details see [0). In this sense the fc^Xcoii(<^) and the transport coefficients T = M,7,C 


stand for the collective response and the coefficients of the Q-mode, respectively, with 
= k‘^T. 

It may turn out that collective motion is so slow, as compared to the dynamics of 
the nucleons, that the transport coefficients can be deduced by expanding the response 
functions around a; = 0. At hrst such a procedure has been applied to the intrinsic 
response (cf. [|^). In this way one may get an approximate solution of the secular equation 
for the position of the poles of the collective response (|3.3|). One may write 


k 


x(^) 


i + x( 0 ))+a,(A 




'1^' 
2 9^2 


= 0 


(3.13) 


uj=0 \ / a;=0 

This form invites to dehne coefficients for friction 7(0) and inertia M(0) by 


7(0) = 


^ doj 




a; 


I a ;=0 


and 


M(0) = 


1 


UJ 


duj 

Id'^x'i.^) 
2 


I a ;=0 


I a ;=0 


(3.14) 


(3.15) 


2 duj'^ ' ^^=0 2 

which in the past have been called the coefficients in ’’zero frequency” limit. The effec¬ 
tive stiffness 0(0) is seen to be identical to the local stiffness of the quasi-static energy 
(mind (^A])), as one would expect to hold true for slow motion, indeed. Incidentally, the 
expression for the inertia can be shown to be a generalization of the one of the cranking 
model to the case of damped motion |^. Unfortunately, for strong damping this M(0) 
becomes very small, and sometimes even negative. For practical applications this does 
not always lead to problems, as usually damping may be so strong that inertia must drop 
out of the macroscopic equations of motion. Nevertheless, this behaviour of the inertia 
M(0) is very unpleasant, but fortunately, one can do better. 

Rather than to concentrate just on the denominator of ( |3.3|) , it is better to take into 
account the full information contained in this expression. As we may recall this was 
done implicitly when constructing the transport coefficients M(a;i),7(a;i) and CiuJi) by 
approximating (p.lO|) by (|3.11|) . An equation like (|3.11|) may be obtained by expanding 


(Xcoii(i^)) ^ in ( p.lOp to second order in cu. In this way one gets 


C 


k‘^Xco\\(u]) 


a ;=0 


x(0) + g(0) 
x(0) 


F(0) 


7 


and 


M 


f ^(Xcoll(tij) ) 

k'^ duj 
1 ^^(Xco11h)~^ I 

2 A)2 5a;2 


-1 


uj=0 


{X{0) + C{0)y 
x^(0) 


■7(0) 


I a;=0 




(3.16) 


(3.17) 


(3.18) 


As compared to the zero frequency limit dehned above, there are two modihcations: 
Whereas all three coefficients obtain a factor of proportionality, only the inertia gets an 
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additional contribution. In most practical applications this proportionality factor is not 
very important, as usually the static response is much larger than the static stiffness: 
x(0) 3>| C(0) |. It is only for small excitations and for deformed oscillator shell models, 
in particular, that the size of | C(0) | becomes comparable to that of x(0). However, 
the additional term in (|3.18|) ensures that the modihed inertia does not drop indehnitely 
anymore with increasing damping. Later on we will demonstrate with the help of numeri¬ 
cal results that ( p.lTD and ( p.l 8 | ) approximate the selfconsistent friction 7 (a;i) and inertia 
Miuji) very well at temperatures T > 2MeV. To distinguish from the zero frequency 
limit (|3.14j) and (|3.15|) we will associate the approximation (|3.17|) and (|3.18|) to ’’the zero 
frequency limit for the collective response function”. 


Finally, we should like to mention that relations similar to the ones given in (|3.16|) - 
( |3.18| ) were obtained earlier in |^, namely for the model case that the collective response 
function just consists of one (approximately Lorentzian) peak. Solving the equations 
for (7(0), 7 ( 0 ) and M(0) given in with respect to C ,7 and M one gets back to the 
equations ( p.l 6 |) -( ^T^ ). 


3.3 Transformation to sharp densities 

In the previous section we found microscopic expressions of transport coefficients for large 
scale motion. For many reasons it is desirable to compare them with those of ’’macroscopic 
models” , like there are the liquid drop model for the inertia and the wall formula for 
dissipation—not to mention stiffness, which we have seen to become identical to the one 
of the static energy anyway, as soon as collective motion becomes sufficiently slow. In |]^ 
many points of principle nature have been clarified about how the macroscopic limit can 
be obtained in microscopic theories, concentrating largely on vibrations around stable 
configurations. Here we like to look at another, more practical, albeit very important 
issue. 

By its very nature, these macroscopic models assume the nuclear density to be constant 
inside some surface at which the density drops from the nuclear matter value down to zero 
within zero range. Commonly this surface is parameterized with the same set of shape 
variables which in our description define the deformation of the mean field, for which in 
the present discussion we have chosen the one collective coordinate Q. Contrary to the 
macroscopic picture, the microscopic calculation leads to a density distribution with a 
soft surface, which in addition depends essentially on the occupation of states. Take the 
simple case of a spherical potential: if a shell with fixed angular momentum I is not filled 
completely the corresponding density distribution will not be spherically symmetric. On 
the other hand, it is clear that the transport coefficients T will depend sensibly on the 
distribution of matter. Thus one needs to employ some specific transformation to relate 
them from one case to the other. This is best done looking at the mean value (or moment) 
of the operator F, as a representative for the average density. 

In the case of our microscopic picture, we have seen the self-consistency condition 
(|3.7|) to imply the relation = k^T between the transport coefficients of the F and 
Q-mode. All we need to do is to search for a similar condition which translates the 
F-motion into that of the sharp density distribution, with the latter being expressed 
through the corresponding shape parameter. Such a relation can be found by applying 
the hypothesis, which we will substantiate below, that the average calculated 

with a density distribution having a sharp surface, can be approximated well by applying 
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an appropriate Strutinsky smoothing to the shell model density. The latter is dehned as 

(F)sharp (3-19) 


where the Uj are smoothed occnpation nnmbers [12|, |13]. How they may be nsed not 


only to calcnlate static expectation valnes like in (|3.19|) , bnt of corresponding response 
fnnctions as well has been stndied in ||2^ and |]T^. The desired relation between {F) 
and Q may now be obtained simply by applying to ( p.l9|) the derivative with respect to 
Q. One hnds 


d{F) 


sharp 


d 


dQ 




I 


jk 3 


^7 ^k 


V., |2 I ^ _ 7 

^ dQdQ ^ 


(3.20) 


which may be nsed to dednce Tgharp = F^^rp- Here, Tsharp represents the transport 
coefficients for the sharp density distribntion, bnt calcnlated for the associated Q-mode. 
Combining this relation with the previons one, we get 


q- 

sharp 


= 


F 

sharp 


= (kx’fTiu,,) 


(3.21) 


As the only one fnrther approximation we have assnmed that both density distribntions 
lead to the same averaged valne of the held operator F. In (p.21|) there appear on the 
very left and on the very right the transport coefficients for the Q-mode, once for the 
sharp snrface of the macroscopic models and once for the collective coordinate specifying 
the mean held. We may add here that the qnantity x"* a physical meaning similar 
to that of a static snsceptibility, hence the choice of this symbol. The only diherence to 
the isothermal snsceptibility is fonnd in using the smoothed occupation numbers fij^Cj) 
instead of nj{ej] T). 

Let us turn now to ’’prove” the hypothesis made. This can be done explicitly for the 
simple case of the deformed oscillator potential. For the more general case we want to ap¬ 
peal to physical intuition, to the extent of accepting the idea that by its very construction 
Strutinsky smoothing commonly does lead to the macroscopic picture. 

Suppose we are given a spheroid whose deformation Q is hxed by the ratio between 
the semi-axes in ^-direction and the one perpendicular to it: Q = zq/hq = ui±/ujz The 
sharp surface S is then given by 


2 ^ 2 ^ 
2/0 4 


(3.22) 


The deformation of the sharp density may be classihed by the average of the quadrupole 
operator Q 2 {f) = 2z^ — calculated as 


r 2 A 

m sharp = J Q2{F)Po{F)dr = —{z^O " Vo) 


(3.23) 


Here, po(h') measures the density of homo- geneous nuclear matter representing A nucleons 
distributed uniformly within this surface. In Fig.7 the quadrupole moment (|3.23|) is 
plotted by the line with stars as function of the deformation parameter Q = lOqjiOz- 
Conversely, the surface (|3.22|) may be interpreted as an equipotential surface of the 
corresponding deformed oscillator potential Q • For such a potential one may compute 

2 For this case the quadrupole operator Qi would be related to our field operator F taken at spherical 
shape by F\q=i = -^mujQQ 2 - 
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single-particle wave functions V5j(r) and from them the moment of the microscopic density 
as 


{Q2)dens = j dr 

3 


(3.24) 


Obviously this moment {Q2)dens depends on the occupation numbers rij. In Fig.7 we show 
curves corresponding to three choices of nj: Dotted line: the occupation numbers are 
hxed at the spherical shape, where the lowest energy states are hlled, and kept constant, 
independent of Q (’’diabatic” situation); dashed curve: at each crossing of states the 
particles and holes are redistributed in such a way that always the lowest states are 
occupied (’’adiabatic” situation); solid curve: here, the smoothed occupation numbers of 
the shell correction method are used. We see that for the diabatic occupation numbers the 
quadrupole moment (|3.24|) differs substantially from that of the sharp density distribution 


{Q 2 )sharp with stars), approximately by a factor of two. Recall that this sharp density 

distribution reflects the deformation of the potential. 

On the contrary, the density computed with the adiabatic occupation numbers follows 
the deformation of the potential on average. It can be said that this feature is one of the 
basic elements of the Copenhagen picture of collective motion 0 (see also [|^). Indeed, 
such a redistribution of particles leads to a consistency between the shapes of the potential 
and the density. For static situations, this leads to the well known relation between the 
occupation numbers in the various directions and the corresponding frequencies of the 
potential. In the treatment of [^] this relation is fulhlled for specihc deformations. In 
Fig.7, the latter correspond to the points where the broken line crosses that with the stars. 
For this broken line it is hardly possible to dehne precisely the derivative of {Q2)dens "with 
respect to the deformation parameter. On the other hand, the condition just mentioned 
can be fulhlled everywhere using the Strutinsky smoothed occupation numbers. In this 
way the derivative is well dehned. Furthermore, the quadrupole moment of the density 
computed with the smoothed occupation numbers practically coincides with that of the 
sharp surface distribution. 


4 Microscopic input 


In this section we are going to specify further details of our treatment of nucleonic dy¬ 
namics. It has already been mentioned that the mere shell model is not sufficient. At 
hnite excitations the effects of collisions cannot be neglected; one even expects them to 
become the more important the higher the nucleonic temperature will be. To treat col¬ 
lisions on the basis of an explicit form of a two body interaction Vresid:i,pi) is hardly 
possible. Therefore we follow another path and parameterize the effect it would have 
on the single particle energies. Details of this method can be found in the publications 
mentioned before, in particular in |T^ (see also |2^, |I8| ). 


4.1 Intrinsic response function 


The Fourier transform of the intrinsic response function given in ( p.4|) can be expressed 
as the sum over single-particle states 


X(w) = 

jk 


(4.1) 
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with 


/ OO , 

n{Q) i^gk{Q)Gj{Q + UJ + ie) + Qj{VL)Qk{S^ - uj - ie)j (4.2) 


Here n{x) is the Fermi function determining the occupation of single-particle levels. The 
Pfc(a;) represents the distribution of single-particle strength over more complicated states. 
It may be parameterized by 


Qk{uj) = 


r(a;) 


(ho; - Cfc - F'(a;))^ 


(4.3) 


in terms of the real and imaginary part of the self-energy S(a;, T) = S'(a;, T) T)/2 

which are assumed to have the following forms: 


n^,T) = 


and 


S'(o;,T) = 


^ {hu - /i)^ + tt^T^ 

ro 1 + ^ 

(ho; — fi) 


—c 


(4.4) 


(4.5) 


2ToVc^ + vr2T2 1 + [{huj - /if + 7r2T2] /c^ 

Both are connected to each other by a Kramers-Kronig relation. The /i in (O)-(O) is the 
chemical potential and the cut-off parameter c accounts for the fact that the imaginary 
part of the self-energy does not increase indehnitely when the excitations get away from 
the Fermi energy. In the present calculation we choose Fq = 33MeV and c = 20MeV. 
The Qk appearing in (|4.2|) is the one-body Green function 


± ie) = 


huj — ek — S'(a;, T) ± iF(a;, T)/2 
which is related to the spectral density Qk by 

Qkiuj) = i{Qk{i^ + i^) - Qk{^^ - i^)) 


(4.6) 


(4.7) 


Details about the evaluation of the integral in (|4.2|) are given in the Appendix. 

For future purpose we want to use this form ( [4.2D of the response function and write 
down a more detailed expression for the friction coefficient in zero frequency limit: 


7(0) = - / 

jk 


(4.8) 


It is obtained (see [|^) by substituting eqs. (|4.1|) - (|4.2|) into (|3.14|) . 


4.2 The problem of ergodicity 

In section 3.1 it was mentioned that the condition of ergodicity (|3.6|) is hard to fulfill in 
the deformed shell model. This statement refers to the study presented in [7^1, where it 
was shown that even collisional damping does not help, at least in the version as used 
to date. One important reason was seen to lie in the fact that our renormalized single 
particle energies have the same degeneracies as those of the pure shell model. But these 
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degeneracies are by far bigger and happen ranch more often than one would expect for 
conhgurations of the compound nucleus. If one believes the latter to be important — 
which should be the case for a hssioning system which has to overcome a large barrier — 
one needs to take special measures to cure the problem just mentioned. 

One way to make the dehciencies apparent is to look at the heat pole which shows up 
either in the correlation function or in the relaxation function <h"(a;). Both are 

related to the dissipative response as 

x"(^) = \ tanh (4,9) 


(see 1^ or []^ for more details as well as for references to the original literature). Both 
functions have a peak at a; = 0, whose width Ft was seen to be twice the single particle 
width calculated at the chemical potential, i.e. Tt = 2r(a; = /i,T). On very general 
grounds, the height of this peak can be seen to be proportional to the difference of 
isothermal and isolated susceptibilities, ~ x(0)- Numerical calculations in |]^ showed 
this height to be large. Rewriting this difference as {x^ — ~x(0)) cind recalling 

that in the nuclear case the difference between the isothermal and adiabatic susceptibility, 
(x"*" ~ X'^'^)) is small [|^], |^, the origin is identihed to come from a large violation of 


ergodicity: x^^ 7^ x(0)- 

This discussion indicates what we can do to cure this problem: Cut the contributions 
of the heat pole to all functions mentioned previously down to the magnitude it would 
have in case the system were ergodic. This means to reduce the height of this peak 
at a; = 0 by the factor (x"*" — X^^)/(x'^'^ ~ x(0))- EH ^ system was studied were 


.ad 


X" ~ X““ vanishes identically, such that the reduction of the heat pole amounted to 
neglect contributions from all states having the same energy. Such a correction can easily 
be done. In expressions like (^4.11) one simply has to restrict the summations in a proper 
way. Here we like to adopt a similar procedure, in the sense of neglecting the influence of 
a hnite difference x"*" ~ At the temperatures considered it will be very small, indeed. 
On the other hand we go one step further and neglect also contributions from neighboring 
states whose energy is hnite but smaller than the collisional width the particles, which as 
we just saw also rehects the width of the heat pole. 

The consequences of such a manipulation are shown in Fig.8 for the intrinsic response 
(left part), the corresponding relaxation function <I)"(a;) (upper right part) and the collec¬ 
tive strength distribution XcoiE*^) (lower right part). All of them have been computed for 
T = 1 MeV. The fully drawn lines correspond to calculations where all matrix elements 
are taken into account. For the dashed curves matrix elements Fjk between states of 
energy difference \{ek — e^)! < F(/i,T) have been discarded. Their contribution to the 
relaxation function is exhibited in the upper right part by the dotted line. This plot 
demonstrates nicely that the heat pole can be associated to a Lorentzian of width Ft 

0 ^ 
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around a; = 0, and thus corresponds to a pole on the imaginary axes 
interesting to note that for the present calculation 99% of this Lorentzian is made out 
only by two pseudo-crossings of single-particle states which takes place close to the Fermi 
energy. By pseudo-crossing we mean a situation where as function of Q two levels come 
close but never cross; one such event is encircled in Fig.3. The big contribution to friction 
in the zero frequency limit of ( |3.14|) (or ( p.lTD ) which results from this heat pole can be 
estimated looking at the upper right part of Fig.8: The slope of the fully drawn line is 
much bigger than the one of the dashed line. On the other hand, contributions to 7(0) 
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from real crossings of levels is very small. This is due to the fact that the matrix elements 
\Fjk\ vanish exactly at the crossing points and are very small in the vicinity of such a 
(real) crossing. 

The influence of the heat pole (in the intrinsic system) on the collective strength 
distribution can be inferred from the lower right part of Fig.8. The fully drawn line shows 
a small peak at very small frequencies, say hoo ~ 0.2 — O.SMefd. Compared to the much 
bigger strength found in the peak at about 1.5 — 2 MeV one is inclined to just ’’forget” 
the small peak when one wants to dehne the transport coefficients. Indeed, it somehow 
looks very natural to associate the bigger peak to the genuine low frequency mode. This 
may be understood as another argument to leave out the contribution of the heat pole 
to the transport coefficients, besides the one involving ergodicity for intrinsic motion. 
All computations to be reported below were done along this line, i.e. contributions to 
response function from the states with the energy difference \{ek — ej)\ < r(/i, T) were 
not taken into account. 

4.3 The influence of collisional damping on nucleonic relaxation 

The whole formulation of our theory is based on the assumption that the nucleonic degrees 
of freedom stay close to a thermal equilibrium. The latter is not fixed, however; rather 
it continuously gets disturbed by collective motion itself, forgetting for the moment a 
possible evaporation of light particles or gammas. It should thus be of interest to have 
some estimate of an appropriate relaxation time. The best candidate for this is offered by 
the ’’generator” of collective motion, namely the F{xi,pi,Qo), which defines the coupling 
of the collective variable to the nucleons. It is predominately this quantity which ’’decides” 
which kind of modes of the intrinsic degrees of freedom get excited. We may recall from the 
discussion in sect.3 the close relation of this F to the nucleonic response function appearing 
in our theory: The x(a;) parameterizes that average ’’excitation” S{F)^ which comes about 
through a change of the collective variable 6Q. If we are just interested to estimate the 
time T after which the S{F)^ has decayed to its static value, we may study the time 
dependent function x{t) given in (|3.4|) . In literal sense it represents the S{F)^ if excited 
by a sharp pulse like SQ(t) oc 6{t). Notice please, that for t > 0 the x(t) is proportional to 
the derivative of the Fourier transformed relaxation function, namely x(t) oc d^"{t) / dt 
(see ( p.4|) and (|4.9|) ). Thus the information contained in x(t) is equaivalent to that of 
<h"(f) up to an additive constant. The latter measures the long time limit of it is 

related to the strength of the heat pole. At the moment we are interested only in the 
behaviour for hnite times. 

In Fig.9. results of two different computations of x{t) are shown. One case just refers 
to the deformed shell model, for the other one collisional damping is taken into account. 
It is clearly seen that only by way of such collisions we may speak of genuine relaxation. 
It is also observed that the latter does not depend much on the shape. However, the 
relaxation time r decreases considerably with increasing temperature. The latter effect 
is expected of course from the very form by which the temperature appears in the single 
particle widths ([4.4| ). From the fully drawn lines of Fig.9 one may deduce for r values 
like 0.5, 0.3, 0.1 fi/MeV for T = 1, 2, 5 MeV. (This may for instance be done by 
approximating the envelope of x(f) by an expontial such that r may be dehned through 
the ’’width” at half maximum.) Later on in sect.5.3 we are going to compare them with 
typical time scales of collective motion, but we may say already here that this microscopic 
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time r is at least one order of magnitude smaller than the one which measures motion in 
the collective variable Q. 

Finally, we should like to mention that the results found from the present computation 


are in accord with those reported in |25 


5 Numerical results for collective transport coeffi¬ 
cients 

In this chapter we will discuss the numerical results for transport coefficients T computed 
along the fission path of As it was already mentioned, for temperatures above 

IMeV the fission path is parameterized in terms of Cassini ovaloids by only one defor¬ 
mation parameter e. Since we prefer to use the parameter ri2 instead of e we have to 
relate the transport coefficients accordingly. This can be done by exploiting the following 
relation, obeying that the a„’s are fixed. 


-‘ri2ri2 — -‘ee 

which simply follows from general properties of coordinate transformations. The deriva¬ 
tive dri 2 /de is obtained by differentiating ( p.4|) , the result is shown by dashed line in 
Fig.2. Recall please, that the deformation dependence of the transport coefficients is de¬ 
fined essentially by the choice of the collective variables. For example, %-^ 2 ri 2 decreases 
with ri2 but is increases as a function of e. Both in the figures as well as in the text 
below we will omit the indices ri2ri2, keeping in mind that the transport coefficients are 
defined with respect to ri2 (even if sometimes they will be shown as function of e). 


dr 


-2 


12 


de 


(5.1) 


5.1 Accuracy of zero frequency limit for the collective response 
function 


The friction coefficient 7(011) and mass parameter M(uji) defined according to (|3.12|) are 
shown in Fig. 10 by solid lines as function of the deformation parameter ri2 for temper¬ 
atures between 1 and 3MeR. They are compared with calculations for which the ap¬ 
proximation (|3.17|)-(|3.18|) is used for friction and inertia, respectively. The latter results 


are marked by dashed curves. As can be seen, this approximation is quite accurate for 
T > 2MeV. This implies that for such temperatures one may avoid the time consuming 
computation of the frequency dependence of collective response function. One may com¬ 
pute friction and inertia directly from (|3.17|) and ( p.l8|) . As for the big fluctuations seen 
at T = 1 MeV we expect them to become much smaller as soon as pairing correlations 
will be taken into account, which shall be subject of future studies. 

For comparison we also show in Fig.10 (by the lines with stars) wall friction 7^.1. 
and the inertia of irrotational flow. According to wall friction 7w.f. is proportional 
to the squared normal velocity m^(s) of the surface, integrated over the nuclear surface. 
Following this may be deduced from the loss of (collective) energy which is given by 


E = -pvp 


ul{s)ds = 7w.f.hi2(t) 


(5.2) 
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where p and vp are the nucleons’ density and Fermi velocity. For axially symmetric shapes 
the surface velocity Un{s) can be expressed in terms of the prohle function p(z, e) from 


eq. (pTTI) as 


Un{z) = Un{z)ri2{t), Un{z) = 


1 dp{z, e) / dru ' 
A(z) de i de 


-1 


A{z) = 


\ 


1 + 


' dp{z,e) 

dz 


(5.3) 

The mass parameter of an incompressible irrotational fluid has been computed as sug¬ 
gested in m- It can be written as 


Mirr = m f ^{s)Un{s)ds 


(5.4) 


with the potential ^{r) for the velocity held expressed by the potential of some ’’surface 
charge” distribution 

1 f u{s') 


e(r) = -i 


ds' 


2ti J \r — f{s' 

The substitution of (p.5|) into the Neumann equation 

= 0, (nVOs = Mn(s) 


(5.5) 


(5.6) 


leads to some integral equation for the density of the ’’surface charge” v[s) which was 
solved iteratively starting with 

(5.7) 


ly- 5 ^ -U„i3 


as zeroth approximation to for details see |^. We have checked that for the par¬ 
ticular case of the shape family (|2.1|) the Werner-Wheeler method turns out to be a 


very accurate approximation to the mass parameter (|5.4|) . Both results coincide within 
the thickness of the lines in the Figure. 


5.2 Temperature and deformation dependence of friction and 
inertia 

In Figs. 11 and 12 the friction coefficient and the mass parameter are presented as function 
of the deformation parameter e for hve different temperatures T = 1 — 5 Meld. They have 
been calculated from the oscillator £t ( p.l2|) , but transformed to sharp densities according 
to (^.21| ). The parameter e is chosen for reasons to be given below. In both hgures results 
of the commonly adopted macroscopic models, namely jw.f. and M^rr are shown by the 
lines with stars. Several observations can be made. 

• To some extent the jw.f, and M^rr can be said to be reached at the higher temperatures. 
The very fact that this statement is more true for the inertia but less so for friction can be 


understood as follows. As shown in |]^ , jw.f. niay be considered the macroscopic limit of 
our model only if such subtleties as collisional damping are left out. Conversely, the high 
temperature limit of the inertia is related to the value of the energy weighted sum, and the 
latter is known to be associated to Mirr- This is true at least when one treats vibrations 
within simple models (see [§> El for the situations of T = 0 and at T 7^ 0). 

• As has been demonstrated in Q , in the general case the evaluation of the sum rule value 
becomes somewhat delicate for collisional damping. The nice feature of reaching the Mj^r 
automatically can be considered a proof of the usefulness of the transformation (|3.21|). 
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• Like in the calculations of [Q, on average friction is seen to increase with T whereas the 
inertia decreases. 

• Again similarly to both transport coefficients decrease with deformation on average. 

• The coordinate dependence shows fluctuations around the average trend, which are more 
pronounced at smaller temperatures. They appear to be even larger than those reported 

0 . 


but they are similar to the results of 


As for the friction coefficient, the 


m 

dependence on deformation is particularly strong around the spherical shape. It reaches 
a kind of local maximum there, or perhaps at a slightly oblate shape. These features can 
be exhibited more clearly for parameterization in terms of e, for which reason we have 
made this choice here. 

Many of these features point to the importance of shell effects, in particular the peak 
for friction around the spherical conhguration, which is clearly visible for T ^ 1 — 2MeV 
but which disappears at T ^ 4 — bMeV. We like to elaborate on this statement by trying 
to split the friction coefficient up into a smooth and a fluctuating part. For this study we 
take the zero frequency limit. 

Let us suppose for a moment that the matrix elements 
of the single-particle energies have some smooth average component JF^(e, e') and the 
oscillating component can be neglected. In this case one may rewrite eq. ([4.8|) in the form 


Fjk\ considered as function 


7(0) ss j 




(5.8) 


The Qki^) in ([4.8| ) are peaked functions with their maximum at hfl = so that the sum 
of Qk over k may be interpreted as the density g{hfl) of single-particle states. As usual 
the latter can be split into the smooth and oscillating components g{e) = g{e) + Sg{e) 
with g{e) = {g{e)) and {6g{e)) = 0 where the brackets stand for an averaging over the 
single-particle spectrum. Inserting this decomposition into (|5.8|) , and noting that an 


integration performed with the bell-like function dn{kl)/dkl can be understood like the 
average introduced above, we will get 


7 ( 0 ) ss hir g) 


(5.9) 


Assuming the oscillating component of the density to be periodic in the energy with 
some period hflo and amplitude Sgo, viz 6g{e) = 6gosm27ie/hflQ, it is easy to convince 
oneselves that {Sg‘^{e)) = {l/2)Sgl. The quantity 6go is determined by the magnitude 


of the shell correction. It is a smooth function of particle number ( see |^) but still 
depends on deformation. It is maximal at that deformation where the shell structure is 
more pronounced. For the Woods-Saxon potential this happens to be so at the spherical 
shape. So it is due to the shell structure that the friction coefficient gains additional 
contributions around the sphere. For the case of the Woods-Saxon potential this specihc 
feature seems to be responsible for the dip one sees in Fig. 10 both for friction as well as 
for the inertia around ri2 ~ 0.5, even at the somewhat larger temperatures of 3 MeV. 
It may be said, however, that such a behavior is not seen in computations performed 
with the model of |^. There the change of the friction coefficient with ri2, for instance, 
resembles more the smooth one given by the wall formula. 


5.3 Collective relaxation times 

It has become customary to parameterize friction in terms of the ratio (3 = y/M, although 
most of the interest has concentrated on friction alone; often the inertia was simply taken 
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to be the reduced mass of the fissioning system. As we have seen above, for a microscopic 
theory both quantities will vary not only with the collective variable but with temperature 
as well. It may be expected, of course, that for the ratio the dependence on shape is much 
weaker than that of the individual quantities. This follows simply from the observation 
that any common, purely geometrical factor will drop out. In Fig.13 we present the 
results for (3 obtained from those for inertia and friction discussed before. Indeed, this 
quantity is essentially constant over the whole deformation region, for all computations 
but T = 1 MeV, a case for which the fluctuations seen in our results are too big because of 
our neglecting pairing. However, there is a marked dependence on excitation: (3 increases 
strongly with T. This is in clear distinction to the result one gets from applying the wall 
formula for friction and that of irrotational flow for the inertia. Interestingly enough, 
these macroscopic estimates lead to some Q-dependence, which in a sense is even bigger 
than suggested by the trend of our results. It can be said that the latter are very close to 
those obtained in the computations with a two-center shell model potential [^]. 

Physically, the inverse of (3 can be interpreted as the relaxation time = M/y) 
to the Maxwell distribution for collective motion. In full glory this feature can only be 
understood looking at the dynamics in collective phase space (see e.g. m)- However, 
one may grasp its content recalling the local equation for average motion of the damped 
oscillator: Mq(t) + 7?(t) -|- Cq(t) = 0. From this equation it also becomes apparent that 
yet another relaxation time can be defined, namely Tcou = '^/\C\. For overdamped motion 
the latter is the only relevant one. Actually, such a situation is given for temperatures 
above 1 — 2 MeV. 

Let us infer this feature comparing the macroscopic relaxation times with the mi¬ 
croscopic r found in sect.4.3. There we got values in the range between r = 0.3 and 
0.1 h/MeV for temperatures between 2 and 5 MeV, discarding the case of T = 1 MeV 
for the moment. Taking for the (3 of Fig.13 typical values like 2 and 8 MeV/fi for T = 2 
and 5 MeV, respectively, we recognize the Tkin to be comparable or only slightly larger 
than T. This means that in this range of temperatures the friction force is so strong that 
it leads to an instant damping or disappearance of the kinetic energy. This is not in con¬ 
tradiction to the basic assumption behind the quasi-static picture. To justify the latter 
all one needs to have is the motion in Q to be slow compared to that for the intrinsic 
degrees of freedom. For this question it is the Tcoii which becomes relevant. Its value can 
be estimated with the help of Figs.6 and 10. Let us concentrate on deformations for which 
the local stiffness does not get zero. Since for T >2 MeV shell effects are not important 
anymore we may estimate C as the one of the liquid drop energy. The right part of Fig.6. 
tells us the \Cldm\ to be of the order of 150 MeV for ri2 larger than about 0.6. In this 
range the value of 7 is about 800 h, as seen from Fig. 10 (for T = 2 MeV ; it is even bigger 
for larger temperatures). For Tcou this implies values of the order of 5 Ti/MeV, which are 
larger than r by more than one order of magnitude. 


6 Summary and conclusions 


In this paper we have applied the single particle model of 0 to describe large scale 
motion at hnite excitations. To this end this model had to be modihed to include effects 
of collisional damping and it had to be adapted to the formulation of collective motion 
in the spirit of the locally harmonic approximation (for a review see |l^). Numerical 


computations have been performed for the transport coefficients of average motion along 
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a fission path. The latter was identihed by the valley in the potential landscape obtained 
for the liquid drop model. It was argued that for the range of temperatures considered in 
the present study this path may be expected to represent fairly well the actual situation. 
First of all, because of the evidence one has from static considerations that Cassini ovaloids 
describe well the shapes of the hssioning nucleus. Secondly, for the large damping one 
expects to be given, the system will be creeping down the collective potential and thus 
will stay close to the line of steepest descent. 

For the transport coefficients values were found which are in accord with previous 
studies, in particular with the ones of 0. This is especially so for their dependence with 
temperature, and to lesser extent for their variation with the nuclear shape. For instance, 
it turns out that typical effects of single particle motion become more apparent here than 
they did in |Q, like there are fluctuations of both friction and inertia with the collective 
variable. To large amount they disappear when building ratios like for the /3 = y/M or 


the rj = 7/ {2^Jm\C\, two quantities which are commonly used to parameterize collective 
dynamics. 

A new development has been achieved with respect to comparisons with macroscopic 
models, like that of irrotational flow for the inertia and that of wall friction. As one knows, 
both are calculated for sharp density distributions which may differ considerably from 
those which correspond to the nucleons’ densities in the shell model. A transformation 
was suggested which allows one to connect both density distributions and thus enables 
one to connect the transport coefficients accordingly. This transformation was found by 
applying a Strutinsky smoothing procedure to evaluate the average static value of the 
relevant one-body operator F{xi,pi,Qo). If is this operator which by a self consistency 
argument is related to the deformation of the mean held. 

Finally, we briehy like to turn to comparisons with experimental hndings. First of all, 
we may mention that the values for (3 shown above concur with the range suggested by 
hssion experiments [^. Two other relevant parameters are the w = ^J\C\/M and the rj 
mentioned above, which for instance appear in Kramers’ famous formula for the decay rate 
of a one-dimensional hssion model: Rk = (a/1 + Ps — hs) ('^m/2vr) exp(—R/T). Here, 
B measures the barrier height, Wm determines the vibrational frequency in the potential 
minimum and rjg is to be evaluated at the saddle. The latter two quantities are shown 
in Figs. 14 and 15, respectively, for both points. They have been calculated by averaging 
the Q—dependent transport coefficients in the neighborhood of the minimum and the 
barrier, where both are those of the T-dependent potential energy. These results agree 
very well with those already shown in ^ . It is remarkable that vj does not change much 
with T, neither at the minimum nor at the barrier. Conversely, like the (3 shown before, 
the Tj too dehnitely increases with temperature. This behavior is in at least qualitative 
agreement with the one found in and the numbers for 7 lie in the range of values 
deduced only recently in from experimental evidence. It is true, of course, that for 
more quantitative analyses one would have to consider both the effects of pairing as well as 
of angular momentum, not to mention the fact that there is some freedom in adjusting 
the two parameters which dehne collisional damping. However, there is probably little 
doubt about the temperature dependence of our coefficients, say above T = 1.5 — 2 MeV. 
This marked change with T is in clear distinction to the macroscopic models mentioned. 
In the light of this feature agreements between experimental hndings and theoretical 
descriptions appear somewhat questionable if they are only based on these macroscopic 
models 0, H]. Likewise the somewhat peculiar behavior of [3 with ri2 suggested in 0 
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is not confirmed by our results (see also 0). It is perhaps fair to say that some of the 
difficulties one still encounters at present when comparing theory with experiment is due 
to the high complexity of the problem itself as well that of the analysis of the data. 
To underline this statement, we just like to take up a point raised in namely that 


common statistical codes evaluate the hssion decay rate not by the transition state result, 
to which Kramers’ formula reduces to for > 0, but by that of the ’’statistical model” 
where vUm is replaced by T. The difference between these two variants can easily be 
inferred from Fig. 14. 
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A Intrinsic response function 


The computation of the intrinsic response function is somewhat involved. Due to the 
frequency dependence both of Fermi distribution n{Lj) and width F(a;) the integral in (4.2) 
can not be calculated analytically. The numerical integration is rather time consuming 
since Pfc(a;) are sharply peaked functions which width varies by two order of magnitude 
depending on the values of e*,. Fortunately the integration in (|4.2|) can be carried out by 
means of residues theorem closing the integration limit in the lower half plane. For this 
one needs to hnd the poles and residues of all the terms in the integrand of (4.2). 

Let us hrst look for the poles of ^^(a; + ie). Substituting (|4.3|), (|T^ ) into ([4.6| ) and 
introducing notation A = \/c^ + = /i ± iA ( [4.6|) can be brought to the form 


— /i ) X 

{(A., - e,){tuv - (A.l) 

It is not difficult to note that the last terms in square brackets cancel with each other 
and both numerator and denominator can be divided by hui — . In this way Qkij^ + 

becomes 


Qki!-^ + ie) = 


hu — /r” 


{huj — huj'^){huj — hcUf.) h^(a;^—o;^) [ a; — o;^ cu — cu 


huj^ — /i ^ /i hujf. 


where hto^ are solutions of the equation 

/ \ _ 

{hu! — huj'j^ ) {hu! — ) = I hu — Ck + i—^ {hu — n ) + 


2Fn 


2FnA 


= 0 


(A.2) 


(A.3) 


namely 




± 





2c^ 


1/2 


(A.4) 


The both poles of ^^(a; + ie) lie in the lower half-plane. The residues of Qkii^ + i^) are 
simple functions of huj^ and as it is seen from (|A.2|) . From the dehnition (|4.6|) it 
is easy to seen that poles and residues of Qkij^ — i^) are complex conjugated to that of 
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^fc(a; + ie). The pole representation for Qk{ui) is easily obtained from ( |A.2| ). Besides the 
poles of Q^{ijj±ie) and Pfe(a;) one should account also for the poles of Fermi function n{ijj) 
in the plane of complex uj (so-called Matzubara frequencies) 

=/i ± i7rT(2?7,-|-1), n = 0,1, 2,... (A.5) 

In principle, the sum extends over inhnite many terms, but in practice the summation in 
(|4.2|) is cut at such frequencies ojn which contribute less than determined by the desired 
accuracy. 
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Figure 2: The relation between ri 2 and e (solid line); the derivative dri 2 /de is shown by 
dashed line. 



Figure 3: The energies of single-particle states (for hxed ; 2 -component of angular mo¬ 
mentum and parity, = 3/2“) as functions of e. The line with stars marks the position 
of the chemical potential computed for T = IMeV. The circle marks the pseudo crossing 
which at e = 0.3 is closest to the chemical potential, see text. 
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Figure 6: The zero frequency limit ( |3.^ ) for the stiffness of the total (left) and liquid drop 
(right) static energies. The dotted, short-dashed, long-dashed and dotted-dashed lines 
correspond to temperatures T = 1,2,3 and AMeV. 



Figure 7: A demonstration of the problem of consistency between the deformation of 
potential and density, exemplified at the deformed oscillator; for details see text. 
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Figure 8: The heat pole contribution to the relaxation and response functions, calculated 
at T = 1 MeV. 
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Figure 9: The time dependent function for the nucleonic response, calculated within the 
pure shell model (dashed curve) and for collisional damping (fully drawn line). 
















Figure 10: The friction coefficient and the mass parameter at hnite frequency (solid lines) 
and for the approximation (|3.17|) - (|3.18|) (dashed lines). 
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Figure 11: The deformation dependence of friction; the dotted, short-dashed, long-dashed, 
dotted-dashed and solid lines correspond to temperatures T = 1, 2, 3,4 and 5MeV. 
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Figure 13: The inverse relaxation time (3 = as function of deformation and 

temperature (indicated in the Figure). 
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Figure 14: The vibrational frequency w = \J\C\/M as function of temperature, calcu¬ 
lated at the potential minimum (marked by circles) and at the saddle point (marked by 
triangles). 



Figure 15: The dimensionless parameter 77 = '-y/{2^\C\M) as function of temperature, 
calculated at the potential minimum (marked by circles) and at the saddle point (marked 
by triangles). 
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